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Abstract 

In the last decades, due to the appearance of many diseases such as SARS and the HlNl flu 
strain, many authors studied the impact of the disease spreading in the evolution of the infected 
individuals using the susceptible-infected-recovered model. However, few authors focused on the 
temporal behavior of the susceptible individuals. Recently it was found that in an epidemic 
spreading, the dynamic of the size of the biggest susceptible cluster can be explained by a 
temporal node void percolation [Valdez et al. PLoS ONE 7, e44188 (2012)]. It was shown that 
the size of the biggest susceptible cluster is the order parameter of this temporal percolation where 
the control parameter can be related to the number of links between susceptible individuals at 
a given time. As a consequence, there is a critical time at which the biggest susceptible cluster 
is destroyed. In this paper, we study the susceptible-infected-recovered model in an adaptive 
network where an intermittent social distancing strategy is applied. In this adaptive model a 
susceptible individual breaks his contact with the infected neighbor with probability a and after 
t\, time units the contact is restored. Using an edge-based compartmental model and percolation 
theory, we obtain the evolution equations of the fraction susceptible individuals in the susceptible 
biggest component. We show that when this strategy is applied, the spreading of the disease 
slows down, protecting the biggest susceptible cluster by increasing the critical time. We also 
study the strategy in a more realistic situation when it is applied after a macroscopic fraction 
of the population is infected. We find that in this case, the strategy can still halt the epidemic 
spreading. Our theoretical results are fully supported by intensive simulations. 

PACS numbers: 64.60.ah, 64.60. aq, 87.19.X- 



2 



I. INTRODUCTION 



Dynamic topologies in complex network models have recently become a subject of 
intensive investigations {l, 2|. While in the past, many dynamical processes were mainly 
developed numerically and analytically on static networks, most of the real networks alter 
their topologies over time. As a consequence, recently many researchers are modeling 
these processes on top of dynamic networks. The topology of the networks can change by 
evolution or by adaptive dynamic l|, l2|. Those networks in which the structure changes 
regardless of the processes taking place on top of them are called evolutive networks. On 
the other hand, networks that alter their topology to mitigate or promote these processes 
are called adaptive networks In adaptive networks, there is a co-evolution 

between the dynamic process and the topology. The study of adaptive networks has 
generated a great interest in many disciplines since there is evidence of many real systems 
with an adaptive topology jli. For example, from the analysis of the gene regulatory 
networks, it was shown that interactions between genes can change in response to diverse 
stimuli, leading to changes in the network's topology jl, 5|. On the other hand, the recent 
pandemics of SARS j^l and HlNl have promoted the modeling of adaptive strategies 
on the contact network to slow down the spread of the epidemics j8-12|. The study of 
these mitigation strategies could provide important information to adopt health policies, 
allowing to characterize the effect of the strategy on the structure of the society and to 
quantify its effectiveness against the spread of an epidemic. 

The most used model that represents these recent diseases is the susceptible-infected- 
recovered (SIR) model 1^ [l^, in which the individuals can be in one of three states, 
susceptible, infected or recovered. In this model, the disease propagates on top of the 
contact network until it reaches the steady state, i.e., when there are no more infected 
individuals. In static networks, the steady state of the SIR model was widely studied 
using percolation tools through a generating function formalism 15l-|l8|. It is known 



that the final size of the disease (fraction of recovered individuals) is governed by a 
control parameter which is the effective probability of infection or transmissibility T of 
the disease. At a critical threshold T = Tc, the disease overcomes a second order phase 
transition with an epidemic phase for T > T^, while for T < the disease consists 
only of outbreaks that reach a small fraction of the population. On the other hand, the 
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adaptive SIR model in the steady state has being explored by few authors. Wang et 



al. 



12| proposed an intervention strategy in which susceptible individuals, induced by 



fear, break the links with their neighbors, with a probability related to the number of 
infected neighbors, regardless of the state of the neighbors. They show that this strategy 
decreases the number of infected individuals and delays the progression of the disease 
compared to the case of non-intervention. Lagorio et al. 19| studied a rewiring strategy 
in which susceptible individuals redirect their links with infected neighbors towards other 
susceptible individuals with probability w. They showed that there is a phase transition 
at a critical rewiring threshold Wc separating an epidemic from a non-epidemic phase 
which can be related to static link percolation. In a recent paper, Valdez et al. ^ 
proposed an adaptive SIR model driven by intermittent connections where the susceptible 
individuals, using local information, break the links with their infected neighbors with 
probability cr for an interval th after which they reestablish the connections with their 
previous contacts. This model, called intermittent social distancing (ISD) strategy focus 
to model the behavior of individuals who preserve their closer contacts during the disease 
spreading. Using the framework of percolation theory, they derived the transmissibility 
To- and found that there exists a critical probability (Tc depending on tb, above which 
the epidemic spread is stopped. They also showed that the ISD strategy, produces a 



■'susceptible herd behavior" below a transmissibility T* |20l . |2l| that protects a large cluster 
of susceptible individuals from being infected. This focus on the susceptible network 
provides a description of the functional network, since the susceptible giant component is 
the one that supports the economy of a society. 



Recently, Valdez et al. 



22| using an edge-based compartmental model (EBCM) 



23 



25] and a generating function formalism, proposed a set of equations which describe the 



evolution of the size of the giant component of the susceptible individuals (GSC) in static 
networks. They found that a time-dependent quantity $5(t), namely, the probability 
that a given neighbor of a root node is susceptible at time t (more details are given 
below), is the control parameter of a temporal node void percolation process involving 
the network composed by susceptible individuals. They showed that the GSC overcomes a 
second order phase transition at a critical time tc above which it is destroyed. This result 
suggests that any mitigation strategy should be implemented before this critical time in 
order to prevent the destruction of the functional GSC. In turn, the critical time can be 
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useful to quantify the effectiveness of the strategy and to compare different mitigation 
measures. 

In this paper, using the EBCM approach and percolation theory, we study the temporal 
evolution of the fraction of the susceptible individuals on adaptive networks following the 
ISD strategy. We find that this strategy protects a giant susceptible cluster by increasing 
the time tc at which the functional network is destroyed. In a more realistic scenario 
where the implementation of the strategy is delayed, we obtain that it also increases 
and find that the dilution of the GSC can also be described theoretically by percolation 
tools. The paper is organized as following: in Sec. [TTl we explain the node void percolation 
that describes the dilution of the susceptible network. In Sec. IIIII we derive the dynamic 
equations of the ISD strategy and in Sec. II VI we present the theoretical and the simulation 
results of the ISD strategy implemented from the beginning of the disease spreading. In 
Sec. |V] we derive the evolution equations when the implementation of the strategy is 
delayed. Finally, in Sec. |Vl] we present our conclusions and outlooks. 



II. NODE VOID PERCOLATION 



In the discrete version of the SIR model, an infected individual infects a susceptible 
neighbor with probability /3 and he recovers after time units since he was infected, where 
tr is called the recovery time. It is known that in this model the control parameter is the 
transmissibility T = 1 — (1 — /?)*'' that governs the final fraction of infected individuals. 
This model can be mapped into link percolation since when the disease traverses a link 
with a transmissibility T, this process is equivalent to occupy that link with the same 
probability T. 

On the other hand, the network composed by susceptible individuals is diluted during 
the epidemic spreading, since when the disease traverses a link, the susceptible individual 
is removed from the susceptible network. However, this process is different than an ordi- 
nary node percolation because the susceptible nodes are not chosen randomly. Instead, 
they are reached by the disease following a link, and therefore, nodes with higher con- 
nectivity are reached with higher probability than nodes with lower connectivity. This 



kind of dilution, that we called node void percolation 



22| leads, in the steady state, to a 



second order phase transition in which the order parameter is the fraction of nodes Si in 
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the GSC, and the control parameter is $5(00), i.e., the probabihty that following a ran- 
dom chosen link, a susceptible node is reached. Moreover, since the susceptible network 
loses the nodes with the highest connectivities, the phase transition at $5(00) = ^sci has 
mean field exponents, as in an intentional attack process independently of the network's 
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topology 

When T increases, 5*1 decreases and therefore $5(00) and T are interrelated. In par- 
ticular, there is an effective transmissibility T* that set the critical value of <^sc at which 
5*1 — 0, and its value can be obtained using a generating function formalism 20, 2l| . 
Denoting the degree distribution as P{k), T* fulfills the equation 



$Sc = Gi[l-r*(l-$5c)], 



(1) 



where T* is the solution of Eq. IQ, ^sc = Gi [{G[y\l)] and = J2k kP{k) / {k)x^-^ 

is the generating function for the excess degree distribution, that is, the degree distribution 
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. In a recent work on static 



of the remaining outgoing links in a branching process 
networks |22j] , it was shown that Eq. ([1]) is valid not only for the SIR model at the steady 
state, but also during the dynamical disease spreading, which imply that the susceptible 
individuals overcome a temporal node void percolation transition. 



III. DYNAMIC FRAMEWORK OF THE ISD STRATEGY 

In the ISD model, initially all the nodes are susceptible except for one node randomly 
infected, that represents the patient zero from which the disease spreads. In this paper, 
we call active links to the links between infected and susceptible individuals. An infected 
individual transmits the disease to a susceptible neighbor with probability /3 and, if he 
fails, the susceptible individual breaks his link with the infected one with probability a 
for tf, time units. After a period both individuals are reconnected and the process is 
repeated until the infected individual recovers at a fixed time t,, > tj, since he was infected. 
If an active link is broken for more than tb = tr — 1 time units, it is restored as a non-active 
link since the infected individual is recovered. This is equivalent to delete that link. In 



this model, the disease spreads with an effective transmissibility T{l3,a,tr,tb) = T„ |20 
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given by, 

T. = ^/3(l-/3r-i(l-ar-i + /3 Yl <Pira,h,a,(3). (2) 

m=l m=ti,+2 

In the first term of Eq. ([2]), — — a)'"~^ is tlie probability tliat an active link is 
lost due to the infection of the susceptible individual at time m given that the active link 
has never been broken in the m — 1 steps since it appears. In the second term of Eq. ([2]), 
P 0(m, tf,, a, P) denotes the probability that an active link is lost due to the infection of 
the susceptible individual at time m given that the link was broken at least once in the 
first m — 1 time units. The probability (f){m, tb, a, which is only valid for m > + 2 is 
given by j20| 



m — 1 

tf,+i 



u=l ^ ^ 

where [• • ■] denotes the integer part function. Thus, the ISD strategy j20| reduces the 
initial transmissibility T towards an effective transmissibility T^r, mitigating the spread of 
the disease. 

Eq. ([2]) can be rewritten as 

T^ = (3f2^m, (4) 

m=l 

with 

= + (5) 

where Qm is the probability that an active link has neither ever transmitted the disease 
nor been broken during the first m — 1 time units, plus the probability that the active 
link was broken at least once in the first m — 1 time units but at time m it is present. 
Here, 0^ = for m < t;, + 2 and is given by Eq. ([3]) for tj, + 2 < m < tr- 

In order to study the evolution of the number of infected and susceptible individuals in 



the ISD model, we use the edge-based compartmental approach |23|-|25] based on a gen- 
erating function formalism. Denoting the fraction of susceptible, infected and recovered 
individuals at time t by S{t), I{t) and R{t), respectively, the EBCM approach lies on 
describing the evolution of the probability that a randomly chosen node is susceptible. In 



order to compute S{t), a link is randomly chosen and then a direction is given, in which 
the node in the target of the arrow is called the root, and the base is its neighbor. Then 
the fraction of susceptible individuals is given by S(t) = "^f^ Pik)6f = Go{6t), where 
Go{x) is the generating function of the degree distribution and 6{t) = 6t is the probability 
that the base does not transmit the disease to the root. In this approach, infection from 



the root to t 
independent 



le base is disallowed, in order to treat the state of the root's neighbors as 
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25j . In order to find 6t we have to compute the probabilities ^R{t) 
and ^i(t) that the base is susceptible, recovered or infected but has not yet transmitted 
the disease to the root. These probabilities satisfy the relation 

+ + <^R{t) = Of (6) 

As $5(t) depends on the second neighbors of the root, we have to disallow infection from 
the first neighbor (base of the root) to the second neighbors, thus 

^s{t) = GM)- (7) 

Note that $s'(t) can be computed as the square root of the fraction of links between 
susceptible individuals, because is the probability that both stubs of a random 

chosen link belong to susceptible nodes. 

Then, using the EBCM approach adapted to discrete times, the evolution of ^^(t) and 
$/(t) in the ISD model is given by deterministic equations, that are only valid above the 
critical transmissibility where there is a macroscopic fraction of infected nodes, 

A^i = -P^i{t), (8) 
A$5(t) = Gi(^i+i)-G'i(et), (9) 

where A is the discrete change of the variables between times t and Eq. (jSj) represents 
the decrease of 9t when the disease traverses an active link. Eq. ^ describes the decrease 
of ^^(t) (see Eq. ([7])). It is important to note that A$5(t) is non-positive, and — A$5'(t) 
represents the new active links. 

The evolution equation of ^lit) is given by, 

A$,(i) = -P^i{t) - A<l>5 + (1 - /3)(1 - a) nt^/^^s{t - tr) + ^(t), (10) 
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where 



^(t) = -a{l-(3)^j{t) + a{l-(3)^j{t-h) 

+ 5^ (7(1 - /3)f]i,._i,+,A$5(t - i - tr), (11) 

with, 

= (1 - P)'^''''"~\l - a)*^-*^+^-i + (12) 

already defined by Eq. ([5]). 

In Eq. flTU]) the first term represents the decrease of $/(t) when the infected base node 
transmits the disease to the root; the second term corresponds to an increasing of $/(t) 
due to the new infections. The third term corresponds to the recovery of the infected 
nodes in active links, which is proportional to: a) the active hnks that appeared tr time 
units ago i.e., A$5'(t — tr) and b) to the probabihty that these active hnks are present 
tr — 1 time units after they appeared and in the last time unit they do not transmit the 
disease neither being broken, i.e., (1 — /3)(1 — cr)^tr- Note that Qtr'^^s{t — tr) is related 
with the active links that appeared tr time units ago, and remain active in their last time 
unit before the nodes recover and the active links dissapear. The last term corresponds 
to breaking and reconnection of active links [see Eq. (ITT]) ]. 

The first term of Eq. ffTTl) . which we call "breaking term", represents the decrease of 
$/(t) due to breaking active links during a period tj,, which is proportional to (t(1 — /3). 
The remaining terms, which we call "reconnection terms", represent the increase of $/(t) 
due to reconnecting active links. While the reconnection term could be thought as minus 
the "breaking term" delayed by tb time units [i.e., a{l — /3)$(t — tb)], that term needs to 
be corrected due to the recovery of some infected individuals that reduce the number of 
active links that could be restored at time t + 1. In Table [J we show a schematic with an 
example of the correction term. Then \&(t) [see Eq. (ITT]) ] represents the net flow between 
the decrease of the probability $/(t) due to breaking active links, and the increase of 
^i{t) due to restoring active links that have being disconnected tb units time ago. 
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TABLE I: Schematic cases of the "reconnection terms" of Eq. f fTT]) for tr = 10 and ti, = 4. 
The first column represents a typical configuration of different cases of the "reconnection 
term". The second column is the probability that an active link cannot be restored at time 
t. In the first column, each cell corresponds to a unit time. The white cells represent the 
time unit at which an active link between the susceptible and the infected node exists, 
the gray ones represent to the disconnection period and the crosses correspond to the 
case where the infected is recovered. The moment at which an active link appears is 
represented by (I), and the moment when this pair, at time t — tf,, breaks for the next 
tf, time units is represented by (II). When some active links break at time t — th at the 
beginning of time t they have to be restored, which is the case of 1). However, other active 
links which break at time t—ti, cannot be restored as active links because the infected node 
recovers. This is shown in the cases 2a) and 2b) where the infected node has recovered at 
i = 2 time units before t (where i is the sub-index of Eq. f[T21) ). The case 2a), corresponds 
to an active link which appears at time t — t^ — i = t — 12 and breaks by the first time 
8 times unit later. The probability for this case is cr(l — /3)*(1 — o")^. Similarly, the case 
2b) corresponds to an active link which appears at time t — tr — i = t— 12 and breaks, 
but not by the first time, 8 time units later. The probability of this case is a{l — /3)08- 

Probability that the active link can- 
Example 

not be restored at time t 

(I) (II) now 

¥^ \ 

1) •• -I I I I I I I I I 

I I I 

cT(l-/3)8(l-a)7 
ail - (3)cPs 



(I) 
t-t-i 



(11) 
t-t. 



now 
t 



2a) 
2b) 



1 

1 




1 

1 
1 
1 
1 
1 


1 ^XIX 



t=10 



|-i=2H 



Note that in Eq. (fTTl) . (1— /3)cr is the probability that an active link breaks. By contrast, 
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^i{t) corresponds to the state of the base node and no information is provided about the 
state of the root node. Therefore, Eqs. flTII]) and (ITT]) correspond to a process wherein the 
infected nodes break their hnks with infected and recovered individuals too. However, the 
breaking of these hnks has no effect on the evolution of I{t) and S{t) because those links 
are not able to transmit the disease. Therefore, Eqs. ffTOj) and f|TT]) are also valid for our 
ISD strategy, where only the susceptible individuals break their links with the infected 
ones. 

The evolution of the fraction of infected individuals is given by 22| 



AI{t) = -AS{t) + AS{t-tr), (13) 

where —AS(t) = — (Go(^t+i) — Goi^t)), as explained above, represents the fraction of new 
infected individuals. The second term represents the decrease of I{t) due to the recovery 
of infected individuals that have been infected tr time units ago. Note that J(t) contains 
the new infected individuals at time t, and also the infected in a previous time. 

The evolution of the fraction of susceptible individuals in the susceptible giant compo- 
nent Si{t) is obtained |22| by solving the equations 

Slit) = Go{9t)-Go{u;t), (14) 
ujt = Ot-GM + Giiut), (15) 

where Ut is the probability that a base individual, which is not connected to the giant 
susceptible cluster, has not yet transmitted the disease to the root at time t; Go{6t) = S{t) 
is the total fraction of susceptible individuals and Go{uJt) is the fraction of individuals 
belonging to finite susceptible clusters. In this paper, we only present our theoretical and 
numerical results for the susceptible nodes, because we are interested here only on the 
functional network. However, we checked that Eq. (1131) is in fully agreement with the 
simulations for the ISD strategy. 



IV. THEORETICAL AND SIMULATION RESULTS 

We apply the ISD strategy on Erdos Renyi networks (ER) with degree distribution 
P{k) = {k)'' exp{—{k))/k\, in which k is the connectivity and (k) is the average degree 
of the network, and Scale-Free networks (SF) with degree distribution P{k) ~ A;"'*' where 
kmin < ^ < kmax and A is a measure of the heterogeneity of the degree distribution. 
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In order to obtain theoretically Si(t), we iterate Eqs. (l8|)- (fT^ with initial condition 
6*0 = 1 — $/(0) = where is the number of nodes in the network. The time 

t is shifted to t = when 1% of the individuals are infected in order to compare the 

n 

simulation and the theoretical curves |23| . 

For the numerical simulations we construct the networks using the configurational 
model 29|, and for the disease propagation we use the algorithm explained in Sec. Illli For 



SF networks, we choose kmin = 2 in order to ensure that the network is fully connected |30| . 
In the simulations, the disease has a transmissibility T without strategy, but it will spread 
with transmissibility [see Eq. ([2])] due to the ISD strategy that is applied since the 
patient zero appears. In our model, at each time step, all the infected individuals try to 
infect the susceptible neighbors through the active links. If they do not succeed, with 
probability cr, the active links will break at the next step. The infected individuals that 
were infected tr time steps ago recover. The updates are done simultaneously after each 
time step. 

In Fig. [T] we plot the time evolution of the size Si{t) of the GSC, for ER and SF 
networks with the same average degree, obtained from the theory and the simulations, 
where we call "size" to the fraction of nodes belonging to a cluster or component. From 
the plots we can see the total agreement between the theory and the simulations. We also 
plot the size of the second biggest susceptible cluster 5*2 (t), in order to show the time at 
which the phase transition develops on the GSC. 



12 




100 200 -25 25 50 75 

t t 



FIG. 1: Time evolution of the size of the susceptible giant cluster Si{t) (black) and the 
second susceptible component 5*2 (t) (red), for = 10^, tr = 20, (3 = 0.05 (T = 0.64) 
without intervention (Q) and with the ISD strategy for cr = 0.50 and tb = I (□, T^- = 0.50) 
and tb = 2 (0, = 0.41) with {k) = 4 in an ER network (T* = 0.46) (a) and a SF with 
A = 2.63, kmin = 2 and k^ax = 500 (T* = 0.40) (b). The symbols correspond to the 
simulations and the solid lines are obtained from Eq. (IT^ . The dashed lines of the 5*2 (t) 
curves are used as guides for the eyes. (Color online). 

As shown in the figure, the strategy increases the critical time tc compared to the case 
without mitigation strategy ((7 = 0), due to the fact that as tf, and a increases, the effective 
transmissibility T^- decreases together with the number of links traversed by the disease, 
protecting the GSC during a greater period of time. In order to quantify tc, we measure 
the position of the peak of 5*2 (t) (see Fig. [1]), which corresponds to the average time at 
which the GSC is destroyed. In Fig. [2] we plot in the plane tb — cr, Ate = tcicr) — tc(o" = 0) 
for different values of cr and tf,, obtained from the simulations. 
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FIG. 2: Phase diagram — a displaying the increasing Ate of the critical time from the 
case without intervention, for different values of a and th, N = 10^, tr = 20, /3 = 0.05 
with (k) = 4 in an ER network (T* = 0.46) (a) and a SF with A = 2.63, fcmm = 2 and 
kmax = 500 (T* = 0.40) (b). The lowest increasing is in black color, and the highest in 
white color. The dashed lines correspond to the case T„ = T* (Color online). 



From Fig. [2] we observe that the finite values of Ate only correspond to > T* 
[see Eq. ([1])]. We can also see that the value of Ate increases when To- — )■ T* from 
above, because below T* the GSC is never destroyed. In SF networks, this increasing is 
less sensitive to the ISD strategy than in ER networks, which is expected since in the 
former, the infected nodes with high connectivities tend to enhance the disease spreading, 
decreasing Si{t) at earlier times. However, there is a large range of values of t^ and a 
that allow to protect the giant susceptible cluster for diseases with not too high values 



of /3 |20|, |3l|. The ability to decrease the transmissibility of the disease spreading is an 
advantage of the ISD strategy for the susceptible network since it gives to the public 
health authorities more time to implement other policies, such as vaccination and to alert 
the health services, before the functional network composed by healthy people disappears. 



V. DELAYED IMPLEMENTATION OF THE ISD STRATEGY 



While most of the adaptive strategies are applied at the beginning of the disease 
spreading, in real situations, a mitigation strategy is rarely implemented when the first 
infected appears 32|, |33| because at this stage the health authorities ignore the severity 



of the disease or they are cautious to declared the alert of an epidemic since any strategy 
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could affect the functional susceptible network. In this section, we study the performance 
of the strategy when it is implemented at a different time after the beginning of the disease 
spreading. 

In order to obtain the evolution equations of the delayed ISD strategy, we consider 
that this strategy is implemented when a fraction x of the population is not susceptible. 

If we denote J{t) = I{t) + R{t), where J{t) is the fraction of the population that was 
reached by the spread (incidence curve), when J{t) < x, the evolution of the disease 
spreading on top of a static network is governed by the Eqs. ([H]), dH]), (fT^ - lH!^ . except for 



^i{t), that is given by [201 



A$,(t) = -(3^jit) - A$5(t) + (1 - T)A^sit - U). (16) 

When J{t) > X at t = r, the ISD strategy is implemented and the evolution of the 
spreading follows Eqs. fl8|)- f|T0|) . f lT3|) - f|T5|) . However, note that when J(r) = x, there 
is a change in the transmissibility, since just after the disease reaches x percent of the 
population, there are still some active links that appeared before t = t, which we call 
"old links". As a consequence, the effective probability that these links would be traversed 
at early stages is not T neither T„. In order to incorporate this issue in the evolution 
equations, first we have to compute = (1 — /3)*''~^A$7(r — + 2;), i.e., the amount of 
old active links which appeared tr — z times units before t = r. This is equivalent to have 
a smaller recovery time z (with 1 < 2; < tr — 1) instead of tr, in the evolution equations 
of the ISD strategy {34 1. As a consequence, the probability that these old links have not 
yet transmitted the disease at time t = t + j with j < z is Qj (see Eq. (^). Here, we 
consider Qj = for j < 1 and j > tr — 1. Then, the equation of $/(t) for the delayed ISD 
strategy needs to take into account the old active links, thus 

+(1 - /3)(1 - a)Qt-r+i'^i,t-r+iil - 61) + ^*{t), (17) 
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where 



-a{l-(3)^j{t) + a{l-(3)^j{t-h) 
+ lJ2^i'^-mu-t,+^A<t>sit-^-tr)] 




1 



(18) 



and Oi and 02 are the Heaviside functions given by 







if t>T + tr-l 



02 



if t> T + U + h-l 



1 — 



1 otherwise 



1 otherwise 



The first two terms of Eq. (|T7|) are the same as those in Eq. ffTOj) . The third term 
corresponds to the recovery of the active hnks which appeared at time t > r — 1. This 
term is different from zero only after tr — 1 time units since r, and has the same form 
as the third term of Eq. (1101) . The fourth term describes the recovery of the infected 
individuals of old active links. Note that this term, has the same form as the previous 
one, in which ^j t-r+i plays the role of A^s{t — tr)- Finally the fifth term corresponds to 
the breaking-reconnection term that is similar to Eq. (fTOl) . 

In Eq. f|T8|) . the first two terms are the same as those in Eq. (fTT]) . The third term 
corresponds to the correction on the reconnection term due to the recovery of infected 
individuals of the active links that appear after t = t. On the other hand, the last term 
corresponds to the correction due to the recovery of the infected individuals of the old 
links. In Fig. [3] we show a schematic of this term. 
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FIG. 3: Schematic of the correction of the restoring term of Eq. (fT8|) at time t, for an 
old active hnk with z = 8 and = 4 that breaks two units time before the infected node 
recovers. Each cell corresponds to a time unit. The white cells represent the time unit at 
which an active link between the susceptible and the infected node exists, the gray ones 
denote the disconnection period and the crosses correspond to the case where the infected 
is recovered. The moment r at which the ISD strategy is implemented, is represented by 
(I), and the moment when this active link, at time t — tb, breaks for the next tb time units 
is denoted by (II). The fraction of active links represented in the plot is proportional to 



In Fig. m we plot Si (t) as a function of t for ER and SF networks, when the ISD 
strategy is applied after the disease reaches a fraction x of the population. 



•t-T-i- 
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FIG. 4: Si{t) for a disease with U = 20, /3 = 0.10 (T = 0.88) and = 10^ when the 
ISD strategy (with a = 0.50 and tb = 10, i.e., T^j = 0.33) is apphed at different moments: 
never (black, Q); after the disease reaches 25% (x = 0.25) of the population (red. A), 
and from the beginning (blue, □) for an ER network with (k) = 4 (a) and a SF network 
with A = 2.63, k^in = 2 and (k) = 4. The symbols correspond to the simulations with 
100 network realizations, and the solid lines correspond to the theoretical solutions. In 
ER network, for the case without strategy Si{t oo) = 0.47. (Color online). 

As shown in the figure, the strategy can protect a GSC, even when it is applied after 
a large fraction of the population is infected. Note that the delayed strategy does not 
change the behavior of the order parameter Si{t) with ^s{t), thus, the critical value of 
at which the giant susceptible cluster is destroyed remains invariant (see Appendix). On 
the other hand, the strategy protects a smaller fraction of individuals in the functional 
network in SF than in an ER networks, because in the former the dilution process of 
the susceptible network is more efficient due to the nodes with high connectivities (see 
Sec. [n]). However, a delayed intermittent connection strategy can still protecting the GSC 
by increasing tc, and thus mitigating the demand of health services and protecting the 
functional network for an extended period of time. 

VI. CONCLUSIONS 

In this paper, we study how individuals based on local information, contribute to halt 
the epidemic spreading implementing an ISD strategy during a disease spreading. This 
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model where the healthy individuals avoid contact with the infected ones intermittently, 
allows to mimic a behavioral response of individuals who try to protect themselves from 
the disease but also try to preserve their closer contacts, such as friendship and working 
partners. 

Using an edge-based compartmental model combined with percolation theory, we show 
that this strategy increases the critical time at which the susceptible giant component is 
destroyed, giving more time to the health authorities to implement other policies against 
the disease spreading. We also study a more realistic scheme in which the strategy is 
delayed, and found that it also protects the functional network. We show that the dilution 
of the GSC for the delayed strategy can also be described by percolation theory. Our 
theoretical framework are fully supported by extensive simulations. 

We believe that to focus on the susceptible network and its topological properties, 
instead that on the infected network, provides a novel description that could be useful 
for the health services to develop new strategies to protect the society and the economy 
of a region. Finally, this complementary view of an epidemic spreading could give a new 
criterion to evaluate the effectiveness of any strategy. 
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Appendix A: Node Void Percolation in a delayed ISD strategy 

A feature of the delayed strategy is that the disease spreads with two different values of 
the transmissibility (T before applying the strategy and T^- during the strategy), and thus 
the final epidemic size cannot be related with a unique value of the control parameter T. 

In order to study the effect of the delayed strategy on the dilution of the susceptible 
network, we analyze the relation between Si{t) and ^s{t) during the disease spreading. 
In Fig. |5]we plot Si(t) as a function of $5(t), obtained from the theory (see Eqs. (l8l)- (fT5l) ) 
and the simulations, when the ISD strategy is applied after the disease reaches a fraction 
X of the population, for ER and SF networks. 
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FIG. 5: Si{t) as a function of ^s{t) for a disease with U = 20, /3 = 0.10 (T = 0.88) and 
= 10^ when the ISD strategy (with a = 0.50 and tt = 10, i.e., T^j = 0.33) is apphed 
at different moments: never (black, Q); after the disease reaches 25% (x = 0.25) of the 
population (red. A), and from the beginning (blue, □) for an ER network with (k) = 4 
(a) and a SF network with A = 2.63, kmin = 2 and (k) = 4. The symbols correspond to 
the simulations with 100 network realizations and the solid black line is the solution of 
Eqs. fl8|)- (fT5|) . The red arrows with label 1 represent the moment at which the delayed 
strategy is applied, and the ones with label 2 correspond to the steady state for the delayed 
ISD strategy. (Color online). 



Even thought there is not a fixed transmissibility, as shown in Fig. [5l the delayed 
strategy does not change the behavior of the order parameter Si(t) with $s(^) found 



in the case without intervention 



22l |. This result is expected since although the speed 



at which the susceptible giant cluster is diluted changes, the size of the GSC at time t 
depends on the amount of susceptible nodes removed by node void percolation, and not 
directly on the transmissibility. As a consequence, the relation between Si{t) and ^s{t) 
holds even for a varying transmissibility. Then, the critical value of at which the 
GSC is destroyed remains invariant. This implies that in a realistic scenario in which 
the transmissibility is varying, if $s(t) is approaching to the GSC is near to being 
destroyed. Thus the value of the distance to the criticallity of the susceptible network 
^s{t) -~^Sc, could be a crucial information for the authorities to decide if more aggressive 
health policies are needed to halt the epidemic spreading and to protect the functional 
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network. 
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